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Abstract 

In this paper, we consider the transmission of confidential information over a K-/i fading channel in 
the presence of an eavesdropper, who also observes K-/i fading. In particular, we obtain novel analytical 
solutions for the probability of strictly positive secrecy capacity (SPSC) and the lower bound of secure 
outage probability (SOP^) for channel coefficients that are positive, real, independent and non-identically 
distributed (i.n.i.d.). We also provide a closed-form expression for the probability of SPSC when the 
p parameter is assumed to only take positive integer values. We then apply the derived results to assess 
the secrecy performance of the system in terms of the average signal-to-noise ratio (SNR) as a function 
of the K and /i fading parameters. We observed that for fixed values of the eavesdropper’s average SNR, 
increases in the average SNR of the main channel produce a higher probability of SPSC and a lower 
secure outage probability (SOP). It was also found that when the main channel experiences a higher 
average SNR than the eavesdropper’s channel, the probability of SPSC improved while the SOP was 
found to decrease with increasing values of k and fx for the legitimate channel. The versatility of the 
K-fi fading model, means that the results presented in this paper can be used to determine the probability 
of SPSC and SOP^ for a large number of other fading scenarios such as Rayleigh, Rice (Nakagami-n), 
Nakagami-m, One-Sided Gaussian and mixtures of these common fading models. Additionally, due to 
the duality of the analysis of secrecy capacity and co-channel interference, the results presented here 
will also have immediate applicability in the analysis of outage probability in wireless systems affected 
by co-channel interference and background noise. To demonstrate the efficacy of the novel formulations 
proposed here, we use the derived equations to provide a useful insight into the probability of SPSC for 
a range of emerging applications such as cellular device-to-device, vehicle-to-vehicle and body centric 
fading channels using data obtained from field measurements. 
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I. Introduction 

With the proliferation of smart deviees, driven by applieations such as the internet of things 
(loT) [[II, device-to-device communications (D2D) [|3 and wearable sensors [|3l-[|5l, privacy 
and security in wireless networking systems have once again been brought to the forefront. The 
wireless medium utilized by each of these applications has an inherent broadcast nature that 
makes it particularly susceptible to eavesdropping. Traditionally, these systems have attained 
secure communications by employing classical cryptographic techniques; e.g., RSA or AES [|6l. 
Unfortunately, these algorithms are entirely disjoint from the physical nature of the wireless 
medium as they assume that the physical layer provides an error-free link. More recently, there 
has been growing interest in information-theoretic security that exploits the random nature of 
the wireless channel to guarantee the confidential transmission of messages ||3. It is widely 
believed that using this type of approach will provide the strictest form of security for physical 
layer communications ||8]| . 

The notion of perfect information-theoretic secrecy, i.e. I{M;C) = 0, where /(■;■) denotes 
mutual information, M is the plane text message and C is its corresponding encryption, was 
first presented by Shannon [|9||. These ideas were later developed by Wyner IITOl . in which 
he introduced the wiretap channel. Under the assumption that the wiretapper’s channel is a 
probabilistically degraded version of the main channel, he studied the trade-off between the 
information rate and the achievable secrecy level for a wiretap channel and showed that it is 
possible to achieve a non-zero secrecy capacity. The secrecy capacity is defined as the largest 
transmission rate from the source to the destination, at which the eavesdropper is unable to obtain 
any information. Csiszar and Komer [fTTlI later extended Wyner’s work to non-degraded channels 
where the source has a common message for both the intended receiver and the eavesdropper 
in addition to the confidential information for the intended recipient. Further developments were 
made in lIT^ . where it was shown that it is possible to achieve secure communication in the 
presence of an eavesdropper over an additive white Gaussian noise (AWGN) channel provided 
the channel capacity of the legitimate user was greater than the eavesdropper’s. The secrecy 
capacity was then shown to be equal to the difference between the two channel capacities. 
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The effect of fading on secrecy capacity was studied in lfT3ll - ll231l . Li et al. [[T3ll . Liang et 
al. ifTdll and Gopal et al. ifTSll characterized the secrecy capacity of ergodic fading channels 
and presented power and rate allocation schemes for secure communication. In [fT^ and [fTTll . 
the secrecy capacity for multiple access and broadcast channels was considered. Barros and 
Rodrigues [fT^ showed that with signal fluctuation due to fading, information-theoretic security 
is achievable even when the eavesdropper’s channel is of better average quality than that of the 
intended recipient. They analyzed the SOP and the outage secrecy capacity for Rayleigh fading 
channels when both the transmitter and the receiver are equipped with a single antenna in the 
presence of a solitary eavesdropping party. A similar analysis for a system consisting of a single 
antenna at the transmitter and multiple antennas at the receiver was presented in ifT^ . This 
was extended to multiple eavesdropping parties in [|^ and over Nakagami-m fading channels 
in [|2TI and [[22]| . It was found that the SOP increases with the number of eavesdroppers and the 
average SNR of the eavesdropper. The outage secrecy capacity was found to increase with the 
Nakagami-m parameter, since an increase in m decreases the severity of fading in the channel. 
More recently, the secrecy characteristics of other commonly encountered fading models such 
as lognormal, Weibull and Rice have also been studied. For example, the probability of SPSC 
of lognormal fading channels was studied in Il2^ . while the probability of SPSC and SOP^ 
of Weibull fading channels was investigated in ^2M . In ll25ll a secrecy capacity analysis over 
Rice/Rice fading channels was conducted and the probability of non-zero secrecy capacity was 
determined. 

While a number of important performance measures for K-/i fading channels [|26]| have pre¬ 
viously been developed such as energy detection based spectrum sensing [|27]| - [|^ and outage 
probability analysis in interference-limited scenarios with restricted values of /i ll29l . to the best 
of the author’s knowledge the secrecy capacity of n-fi fading channels has yet to be reported in 
the open literature. Due to the equivalency pointed out in [|30l, the results presented here will also 
have immediate applicability in the analysis of outage probability in cellular systems affected 
by co-channel interference and background noise, and the calculation of outage probability in 
interference-limited scenarios. Indeed the new equations proposed here provide an alternative. 
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more general result than those presented in If29]| (wherein the outage probability in interferenee- 
limited seenarios is restrieted to partieular values of the /r parameter) and allow the ealeulation 
of the relevant eapaeity formulations for arbitrary, real, positive values of the /i parameter. 
Motivated by this, we analyze the seereey eapaeity of K-/r fading ehannels in whieh we assume 
the eavesdropper to be passive and the ehannel state information (CSI) of the eavesdropper and 
the intended reeipient are not available at the transmitter. 

The main eontributions of this paper are now listed as follows. Firstly, we derive novel 
analytieal expressions for the probability of SPSC and SOP^ over K-/i fading ehannels for 
eoeffieients that are positive, real, i.n.i.d. for both the legitimate and non-legitimate parties. 
Seeondly, we provide an exaet elosed-form solution for the probability of SPSC over K-/i fading 
ehannels when the /i parameter is assumed to take positive integer values. These expressions have 
been subsequently verified by reduetion to known speeial oases. Thirdly and most importantly, 
beoause the fading model ll26ll oontains a number of other well-known fading models 
as speeial oases, the novel formulations presented in this paper unify the seereey eapaeity 
of Rayleigh, Rioe (Nakagami-n), Nakagami-m and One-Sided Gaussian fading ehannels, and 
their mixtures. Therefore they oan be used to provide a useful insight into the seereey eapaeity 
of eavesdropping seenarios whieh undergo generalized fading oonditions. Finally, we provide 
an important applioation of these results to estimate the probability of SPSC of a number of 
emerging wireless applieations sueh as eellular devioe-to-deviee, vehiele-to-vehiole and body 
eentrio oommunieations using data obtained from real ehannel measurements. 

The remainder of this paper is organized as follows. Seetion II provides a brief overview of the 
K-yii fading model. Seetion III explains the system model while Seetion IV provides the derivation 
of novel analytieal and elosed form expressions for the probability of SPSC and SOP^. Seetion 
V provides an overview of the model parameters required to obtain the seereey eapaeity of the 
eommon fading models derived from the K-/i fading model; this is followed by some numerieal 
results. Seetion VI diseusses some of the applieations of this paper. Lastly, Seetion VII finishes 
the paper with some eoneluding remarks. 
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II. An Overview oe the K-/i Fading Model 


The K-jji fading model was originally conceived for modeling the small-scale variations of 
a fading signal under line-of-sight (LOS) conditions in non-homogeneous environments |[2^ . 
The K-jjL fading signal is a composition of clusters of multipath waves with scattered waves of 
identical power with a dominant component of arbitrary power found within each cluster. The 
received signal envelope, R, of the fading model may be expressed in terms of the in-phase 
and quadrature components of the fading signal such that eq. 6 ] 

= {X, + (11 + qif ( 1 ) 

i=\ i=l 

where p is the number of multipath clusters, Xi and 11 are mutually independent Gaussian 
random processes with mean E[Xi] = £'[11] = 0 and variance E [Xj^] = £ [TJ 2] = (i.e., 

the power of the scattered waves in each clusters). Here pi and qi are the mean values of 
the in-phase and quadrature phase components of multipath cluster i and d? = 

Letting 7 represent the instantaneous signal-to-noise-ratio (SNR) of a K-p fading channel, then 
its probability density function (PDF), ( 7 ), is obtained from the envelope PDF given in If2^ 
eq. 11 ] via a transformation of variables as. 


M+l n-1 -m(1 + w)7 

. / N _ /^(l + «) " 7 2 e 7 
E/ (7) —-M-i _M+i- 


2 'y 2 



K (1 -1 k ) 7 
7 


( 2 ) 


where k > 0 is the ratio of the total power of the dominant components {(R) to that of the 
scattered waves {2pa'^) in each of the clusters, /i > 0 is related to the number of multipath 
clusters and is given hy p = where E(-) and V(-) denote the expectation and variance 

operators, respectively, 7 = E ( 7 ), is the average SNR and /„(•) is the modified Bessel function 
of the first kind and order n. As indicated in [|26l, it should be noted that the p parameter can 
take non-integer values which can be due to the non-zero correlation between the in-phase and 
quadrature components of each cluster, non-zero correlation between the multipath clusters or 
non-Gaussian nature of the in-phase and quadrature components. The cumulative distribution 
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function (CDF) of 7 can be obtained from [|2^ eq. 3] as, 



(3) 


where Q. (•, •) is the generalized Mareum Q-funetion defined in lISTl eq. 4.60] as, 



(4) 


The K-jji distribution is a generalized fading model whieh also eontains as speeial oases other 
important distributions sueh as the Rioe (/r = 1 and n = K), Nakagami-m (k —)■ 0 and /i = m), 
Rayleigh (/r = 1 and k ^ 0) and One-Sided Gaussian (/i = 0.5 and k —)■ 0). 


III. The System Model 


Consider the system model of seeure data transmission shown in Fig. [T] The legitimate 
transmitter, Alioe (node A), wishes to oommunioate seoretly with the legitimate reoeiver. Bob 
(node B), while a third party. Eve (node E), is attempting to eavesdrop. We assume that the 
main and eavesdropper’s ohannels both experienoe K-jj, fading. Alioe wishes to send a message, 
= [te(l),te(2),t(;(3).. .w{k)] to Bob. At the transmitter, the message is eneoded into a 
eodeword, = [a;(l), a;(2), a;(3)... x{n)\, for transmission over the ohannel. The signal reeeived 
by Bob and Eve oan be written as. 


Vm ii) = hu ii) X [i) + hm ii) 


(5) 


Ve (*) = hE (i) X (i) + He (i) 


( 6 ) 


where and hE(i) are the quasi-statio K-/i fading ooeffioients of the main and the eaves¬ 


dropper’s ohannels respeotively (i.e., hM(i) = hu'^i and /^^(i) = and nui}) and are 


the zero-mean oiroularly symmetrio oomplex Gaussian noise random variables with unit varianoe 
at Bob and Eve respeotively. 


We assume that the fading ooeffioients of Bob and Eve’s ohannels, although random, are 
eonstant during the transmission of an entire eodeword and independent of eaeh other. Letting 
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Alice - Node A 


Bob - Node B 



V_ 


Eve - Node E 


Fig. 1. The proposed system model. 


P, Nm and Ne represent the average transmit power, noise power in the main ehannel, and noise 
power in the eavesdropper’s ehannel respeetively, then, the eorresponding instantaneous SNR’s 
at the legitimate receiver and the eavesdropper are given by, 7 m = and 7 e = 


Ne 


while the average SNR’s are given by, 7 m = ^ and 7 ^ = K Now let us consider 

the channel components of Bob and Eve which are assumed to be i.n.i.d. random variables with 
parameters {km, 7 m} and {ke, /tb, 7 b}, respectively. The PDF’s of 7 m and 7 ^ can be 
re-written from (l2l) as, 

MM + l MM-l 

„ I N /^m(1 + ^ 7m 2 e I nm{^ + nm)'Jm , 

JiM 17m j --177+1-1 ^dM\l -=- I (/) 


_ MM + 1 

Nm 2 7 m 2 


7m 


f \_dE0- + NE) 2 7b 2 e /kb(1 + kb)7b , 

JiE \1 e) - - 7711 -77+1-Vb-1 1 ^dE\l -=- I • (o) 


Ne ^ 7 e ^ e^^Ei^E 


lE 


Likewise, the CDF’s of 7 m and 7 b can be re-written from ([3]) as. 


-^7m (7m) — 1 Q 


I^M 


F'te (7e) - 1 - Q 




/-- 2{1 + Km) IJ-mIm 

y ^nm^m, \ -^- 

V 7m 

/2a+K|W 


(9) 


( 10 ) 
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IV. Secrecy Capacity in k-h Fading Channels 


In this section, we derive an analytical expression for the SOP^. We also arrive at both an 
analytical and closed form expression for the probability of SPSC over fading channels. 
Herein, it should be noted that we adopt the notation M/E to describe the fading conditions 
experienced by the main channel (M), i.e. the intended recipient, and the eavesdropper (E). 
For example k-jj, / k-jj, indicates that the main and eavesdropper’s channels are both subject 
to i.n.i.d. K-jj, fading. Of course, because of the generality of the K-jj, fading model, M and E 
can be readily interchanged with the Rayleigh, Nakagami-m, Rice (Nakagami-n) and One-sided 
Gaussian models which all appear as special cases of this fading model. 

In passive eavesdropping scenarios, where the CSI of the eavesdropper and the intended 
recipient are not available at the transmitter, perfect secrecy is not guaranteed. Hence, we adopt 
the probability of non-zero secrecy capacity and SOP as useful secrecy metrics to characterize 
the system performance. For the system under consideration, the capacity of the main channel 
is given by Cm = log2 (1 -f 7 m) and the capacity of the eavesdropper’s channel is given by 
Ce = log 2 (1 + 7e)- From flSll, we define the secrecy capacity, Cg, for one realization of the 
SNR pair ( 7 ^, 7 e) of the quasi-static complex fading wiretap channel as. 



log 2 (1 4- 7m) - log 2 (1 -f 7i?), (7 m > 1 e ) 
3 , (7m < 1 e ) 


( 11 ) 


A. SOP Analysis 

The outage probability of secrecy capacity is the probability that the instantaneous secrecy 
capacity falls below a target secrecy rate Rg {Rg > 0 ) and is defined in ^ as. 


Rout (Rg) = Pr {Cg < Rg) . 


( 12 ) 


Performing the necessary mathematical manipulations, we obtain. 


Pout (Rs) — Pr [7M < (1 + 7e) ~ 1 ] 


(13) 


which can then be expressed as, 


p (p \ + K,e) ^ 

out \ S) 1-^E — ^ 


M£; + l 


E~^ _ ^E'^ 

ke ^ jE ^ 


tip,-! -M£(1+ks)7£ 

lE~^e ^ ( 2he\ 


Ke (1 + E,e) Je 
lE 


X 


1 - Q,„ (^s/2^, j j djE. (14) 


Proof: See Appendix |S ■ 

At present, due to the eomplieated form of the integral eontained in (fT4l) it is not possible to 
obtain a elosed-form expression for the SOP, therefore, we derive the lower bound of SOP as 
follows ll2^ . 


Pout [Rs) = Pr [iM < (1 + 7 e) - l] 

> SOP^ = Pr [-iM < e^^-iE] . 

Consider the Mareum-Q funetion in (fT4l) . This ean be re-written as. 

Now, using dH]), dH) and dUl) the lower bound of SOP is obtained as. 


(15) 


(16) 


MS + l 00 


SOP^ = 


L _ /^e(1 + ^e) ^ 

—1 _ A^E + 1 

Ke -fE ^ et^Ei^E 


7E^~e ^ I 2fiE] 


Ke (1 + Ke) -fE 

-fE 


X 




The solution for (fTTI) can be obtained via the following proposition. 


Proposition 1. For the arbitrary real and positive pm, Fe, km and ke (fTTI) can be evaluated 
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as, 


SOP^ = 


Mg + l 




{V OCePe^ 


/ifi —l+2fc 


ttE 2 k=o 


T{fiE + k)r 1 

\Pe) 


Mg + l 

/3e ^ 


aE 2 


X 


OO CXD 


EE 


k\T {fiE + k) 

{V O'ePe)'^^ ^ {‘^O/mY T {flE + k + fiM + ^) 


^ k\r {fiE + k) /!r (/iM + 1) 2 'e"A^ + k) (/^^e^s + 

/3e 




X 2-Pi ( 1, /iij + A; + fiM -\- k, i~ie k -\-1] 


13 + /3e 


(18) 


where a — 1/7m, b — I/je, 13m — {K,M + l)anM, (3 e — {nE + l)bnE, otM — i^Mh-M, c^e — E,Ej^E 
and 2 -Pi(- , • j ■ j ■) is the Gauss hypergeometric function 


Proof: See Appendix |Bl 


B. SPSC Analysis 

In this subseetion we examine the eondition for the existenee of strietly positive seereey 
eapaeity. This oeeurs as a speeial ease of the seereey outage probability when the target seereey 
rate, Rs = 0. Aeeording to |[8l the probability of non zero seereey eapaeity is defined as, 

Pq = Pr (Cs > 0) = Pr (7m > 7e) (19) 


Performing the neeessary mathematieal manipulations, the probability of SPSC ean be expressed 
as. 


Mg+i 


P. = 


/i£;(l + Ke) ^ 
Ke 2 je 2 


yE^e ^ ( 2pE\ 


Ke (1 + Ke) 7£ 
lE 


X 




( 20 ) 


Proof: See Appendix O ■ 

As ean be seen, the solution obtained for (l20l) depends on the parameters pm and pe and 


thus the following propositions are valid. 
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Proposition 2. For the arbitrary real and positive pm pe, (I201) can be evaluated as, 


f^E 2 (\/ ocePe)^^ (pe + k + Pm + 1 ) 

aE^e-^ ^ klT {pE + k) l\T {pM + 1) 2'e“M + k) (/3m + 


X 2 F 1 


^1, pe + fc + Pm F 1] pE + fc + 1; 


Pm + Pe ) 


( 21 ) 


where a = I/^m, b = I/^e, Pm = («M + l)a/^M, Pe = {nE + l)hpE, aM = i^mFm, ^e = f^EdE 
and 2-^1 (■ , • j ■ j ■) is the Gauss hypergeometric function /f32l/ . 


Proof: See Appendix |Dl 


Proposition 3. For integer values of pm and p e, an exact closed form solution for (l20l) is 
obtained as 


Pq = 1— P + exp 


2R 




m=—fi 

\—^ I V k 
X < ^ I 

k=i \ k + m 




‘-E 

i=i 




m 


( 22 ) 


where, 


P= Qi 


2KEapEdM (1 + ^m) 


2KMbpEdM (1 + K,e) 


) 


{b (1 + He) Pe) + (a (1 + km) Fm) V ((^ (1 + ^e) Fe) + (a (1 + km) Fm) ) 


6 (1 + ke) Fe 


exp 


X Jr 


{b (1 + Ke) Fe) + (a (1 + km) Fm) 

f 2 peFms/ abKM^E (1 + ^ e ) (1 + k,m) 

I (& (1 + «^e) Fe) + (a (1 + Km) Fm) 


OKeFmFE (1 + ^ m ) + bKMFMFE (1 + ^ e ) 
{b (1 + Ke) Fe) + (a (1 + km) Fm) 


r = A = ^2keFe] B = y/2KMFM] F = Fe - P, v = pm - 1] and R = r + r ^ 


Proof: See Appendix El ■ 

V. Special Cases and Numerical Results 

In this seetion, we verify the novel analytieal and closed-form expressions of SPSC derived 
above by reducing them to a number of known special cases. We also use the formulations to 
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provide a useful insight into the behavior of the SOP^ and SPSC as a funetion of the fading 
parameters of the legitimate and eavesdroppers ehannels. 


A. Some Special Cases 

As diseussed previously, beeause of the generality of the K-p fading model the results presented 
here eneompass the probability of SPSC for a wide range of fading channels. For the reader’s 
convenience a detailed list, along with the corresponding parameter values, is presented in Table 
m We now discuss a few examples from Table HI 

1) Rice / Rice and Rayleigh / Rayleigh: As shown in Table HI to obtain the probability of 
SPSC for the case when both the main channel and eavesdropper’s channel undergo Rician 
fading (i.e. a Rice / Rice fading scenario), we substitute pm = Pe = ^ into (l2TI) and/or (l22l) in 
which case these reduce to 


(1 + ke) 

rn — - 


oo oo 


{I + Ke) 

'i>—n i—n 


f^E (1 + i^e) b 


2k 


y (2KAf)^ 

— — k\T{l + k) /!r(l + /) 2 'e^M 


X 


^ k=0 1=0 * V ' / • \ / 

- ((>=« + 1) (2 + + 0 ^pfj^_2 + k + l-,2 + k,- -(^£±11^- 

(1 +/c) ((km + 1) a + (ke + 1) 6) \ {km + ^) a +{ke + ^)bj 

(2 


and 


Po = I-Q 



2ke(i (1 + Km) 


2km^ (1 + ^i?) 


6 (1 + ke) + o (1 + Km) V ^ (1 T ^e) + u (1 + Km) 


fe (1 + Ke) 


6 (1 + Ke) + fl (1 + Km) 




2 a / abK,MK,E (1 + i^e) (1 + ^ m ) 
6 (1 + Ke) + fl (1 + Km) 


(24) 


which are in exact agreement with the result reported in ll25l eq. 10] and are illustrated visually 
in Fig. [2l Of course letting km = ke = 0, then 


Po = 


7 m 

7 m + 7 b 


(25) 


which matches exactly with that given in ifT^ eq. 5] for a Rayleigh / Rayleigh fading scenario. 
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TABLE I 

Parameter Substitution to Obtain Secrecy Capacity of Special Cases Derived from the k-h Fading Model 


Seenario 

Parameters to be substituted in (|2TI) and/or (l22l) 

Riee / Riee 

k'M — k'E — 1 j ^ 0 and > 0 

Nakagami-m / Nakagami-m 

/^M — k^E ^ ^ \ I^M —^ 0 t^E —^ 0 

Rayleigh / Rayleigh 

/iM = /is = 1 ; ^ 0 and ke —)■ 0 

One-Sided Gaussian / One-Sided Gaussian 

Hm — k'E ~ 0.5 ; Kjm —^ 0 and ke — )■ 0 

K-yii / Riee 

/lat ^ 0 ; Kjm ^ 0 ; yiiE = 1 and ke ^ 0 

K-fi / Nakagami-m 

fiM ^ 0 j > 0 ; /iE > 0 and ke —^ 0 

K-gi ! Rayleigh 

Hm > 0 ; Km > 0 ; /iE = 1 and ke —)■ 0 

K-fi / One-Sided Gaussian 

jj^M ^ 0 ; Km ^ fiE ~ 0.5 and ke —^ 0 

Riee / 

/lm = 1 ; > 0 ; /le > 0 and ke > 0 

Riee / Rayleigh 

k'M — 1 j > 0 ; yU-E = 1 and ke —^ 0 

Riee / Nakagami-m 

fiM = 1 ; i^M > 0 ; /le > 0 and ke —)■ 0 

Riee/ One-Sided Gaussian 

k'M ~ 1 5 ^ 0 ; fiE ~ 0.5 and ke —^ 0 

Nakagami-m /«:-/! 

fiM > 0 and km —)■ 0 ; /ie > 0 and ke > 0 

Nakagami-m / Riee 

fiM > 0 and Km —)■ 0 ; /ie = 1 and ke > 0 

Nakagami-m / Rayleigh 

fiM > 0 and Km —^ 0 ; /ie = 1 and ke ^ 0 

Nakagami-m / One-Sided Gaussian 

fiM > 0 and Km —> 0 ; /le = 0.5 and ke 0 

Rayleigh / n-jj, 

k'M — 1 j ^ fiE ^ 0 and Ke ^ 0 

Rayleigh / Riee 

fiM = 1 ; k,m fiE = 1 and Ke > 0 

Rayleigh / Nakagami-m 

k'M ~ 1 5 ^ 0 ; yU-E ^ 0 and ke —^ 0 

Rayleigh / One-Sided Gaussian 

k'M ~ 1 j ^ 0 ; fiE — 0.5 and ke —^ 0 

One-Sided Gaussian / n-fi 

fiM — 0.5 ; Km — y 0 ; yiiE ^ 0 and ke ^ 0 

One-Sided Gaussian / Riee 

fiM — 0.5 ; Km —^ 0 ; /le ~ 1 and ke > 0 

One-Sided Gaussian / Rayleigh 

fiM — 0.5 ; Km —^ 0 ; yiiE = 1 and ke —^ 0 

One-Sided Gaussian / Nakagami-m 

fiM = 0.5 ; Km —)■ 0 ; yiiE > 0 and ke ^ 0 


2) Nakagami-m / Nakagami-m: By letting km —^ 0 and ke ^ 0 into (l2TI) and/or (l22l) . we 
obtain the probability of SPSC for the seenario where both the legitimate and non-legitimate 
users ehannels undergo Nakagami-m fading. As shown in Fig. [2l our results have been eompared 
with that reported in [1211 eq. 8]. It should be noted that the expression proposed in eq. 8] 
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Fig. 2. The probability of non-zero secrecy capacity versus 7 m. Lines represent equation (ED and (|22j for the special case of 
Rice (fiM = = 1; Km = 15, ke = 12); Nakagami-m (km —>■ 0 and ke —>■ 0; /tm = ^^E = 2) and Rayleigh fading {km —>■ 

0 and acb —>■ 0; = fJ-E = 1). Triangle markers represent EH eq. 10] with Ka = 15, Kb = 12 for Rice; circle markers ED eq. 

8] with m = 2; N = 1 for Nakagami-m and square markers (H eq. 5] for Rayleigh fading. 


is valid only for identical fading parameters of the main and the eavesdroppers channels and 
for integer values of the shape parameter, m (or equivalently /i when k —)■ 0) when a single 
eavesdropper is considered whereas the equation proposed here is valid for any positive real 
value of the /i parameter. As shown in Fig. |2l for the case of integer and /xg, the results 
presented here are in exact agreement with those presented in ll2T]| . 

3) Other Fading Scenarios: In a similar manner, the probability of non-zero secrecy capacity 
for several different fading combinations, most of which have not been reported previously in 
the open literature, can be obtained from Table HI Fig. [3] shows the shape of probability of SPSC 
for a selection of these scenarios, namely the Nakagami-m / One-Sided Gaussian {km —^ 0 ; 
Hm = 1.2 ; /i£; = 0.5 and ke ^ 0), k-^ / One-Sided Gaussian (hm = 1-2 ; km = 5 ; /i^ = 0.5 
and Ke ^ 0) , One-Sided Gaussian / One-Sided Gaussian (/xm = /xb = 0.5 ; km 0 and ke 
—)■ 0) and Rice / Nakagami-m (/xm = 1 ; i^m = 5 ; /xe = 1.2 and ke ^ 0 ). 

B. Numerical Results 

Here, we discuss the behavior of Pq and SOP^ as a function of parameters {km, Fm, 7m} 
and {ke, Fe, 7e}- Let B represent the ratio of the average SNR of the main channel to that of 
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7m (dB) 

Fig. 3. The probability of non-zero secrecy capacity versus 7 m for the special case of Nakagami-m / One-Sided Gaussian (km 
0; fiM = 1.2; fj,E = 0.5 and ke ^ 0), n-fi / One-Sided Gaussian (km = 5; fiM = 1.2; ke ^ 0 and ^e = 0.5), One-Sided 
Gaussian / One-Sided Gaussian (fiM = fJ-E = 0.5; km — >■ 0 and ke —>■ 0) and Rice / Nakagami-m (km = 5; fiM = 1; k_b —>■ 
0 and ^E = 1.2) 


the eavesdropper’s channel i.e, B = ^m/ie and r represent the ratio, ke/k,m- For the figures 
shown in this section, we set Rs = I dB, /im = 1-4 and = 1-2. Four example profiles for 
the probability of SPSC and SOP^ are illustrated in Fig. 01 


From these figures, we observe that the probability of SPSC increases and the secrecy outage 
probability decreases as B increases. This is because larger values of B indicate that the signal 
quality of the main channel is better than that of the eavesdropper’s channel. From Figs. Ufa) 
and (b), we observe that the probability of SPSC is non-zero even when B < 1. Furthermore, 
Pq increases as km, the main channel’s fading parameter, increases, while for a fixed km with 
7m < lE, Pq increases as ke decreases. It can be seen from Figs. Htc) and (d) that the SOP 
decreases as km increases. Additionally, for a fixed km and when B < 1, the SOP decreases 
as Ke decreases. Similar observations are made when km and ke are set and fiM and /i^; are 
varied. 
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Fig. 4. The probability of non-zero secrecy capacity with (a) km = 4 and (b) km = 10, and the secrecy outage probability 
with (c) Km = 4 and (d) k,m = 10. 


VI. Applications of SPSC for Devices Operating in Fading Channels 

To illustrate the utility of the new equations proposed here, we now analyze the probability 
of SPSC for a number of emerging applieations such as cellular device-to-device, body area 
network (BAN) and vehicle-to-vehicle (V2V) communications using channel data obtained from 
field trials. For all of the measurements conducted in this study, we considered a three node 
system which consisted of Alice, herein denoted node A, which acted as the transmitter and 
also Bob (node B) and Eve (node E) which acted as the receivers. Each of the nodes. A, B 
and E, consisted of an ME5805 transceiver, manufactured by REMD. The transceiver boards 
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were interfaced with a PIC32MX which acted as a baseband controller and allowed the analog 
received signal strength (RSS) to be sampled with a 10-bit quantization depth. For all of the 
experiments conducted here, node A was configured to output a continuous wave signal with a 
power level of -i-21 dBm at 5.8 GHz while nodes B and E sampled the channel at a rate of 1 kHz. 
The antennas used by the transmitter and the receivers were -1-2.3 dBi sleeve dipole antennas 
(Mobile Mark model PSKN3-24/55S). 

A. Device-to-Device Scenario 

The first set of measurements considered cellular device-to-device communications channels 
operating at 5.8 GHz in an indoor environment. The experiments were conducted in a large 
seminar room located on the first floor of the ECIT building at Queen’s University Belfast in the 
United Kingdom. The seminar room had dimensions of 7.92 m x 12.58 m x 2.75 m and contained 
a number of chairs, some desks constructed from medium density fiberboard, a projector and 
a white board. For these measurements, the antennas were housed in a compact acrylonitrile 
butadiene styrene (ABS) enclosure (107 x 55 x 20 mm). This setup was representative of the 
form factor of a smart phone which allowed the user to hold the device as they normally would 
to make a voice call. Each antenna was securely fixed to the inside of the enclosure using a 
small strip of Velcro®. The antennas were connected using low-loss coaxial cables to nodes A, 
B and E. 

The experiment was performed when the room was unoccupied except for the test subjects 
holding the devices. As shown in Eig. [5ta), this particular scenario considered three persons 
carrying nodes A, B and E who were positioned at points X, Y and Z respectively. The three 
test subjects using nodes A, B and E were adult males of height 1.72 m, 1.84 m and 1.83 m; 
mass 80 kg, 92 kg and 74 kg, respectively. During the measurement trial, all three persons were 
initially stationary and had the hypothetical user equipment (UE) positioned at their heads. The 
persons at points Y and Z were then instructed to move around randomly within a circle of 
radius 0.5 m from their starting points while imitating a voice call. A total of 74763 samples of 
the received signal power were obtained and used for parameter estimation. 
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Figs. a) and (b) show the empirical PDF of the signal envelope for Bob and Eve compared to 
the n-fi PDF given in [l2^ eq. 11] for the D2D channel measurements. All parameter estimates for 
the K-fx fading model were obtained using the Isqnonlin function available in the Optimization 
toolbox of Matlab along with the n-fx PDF given in [|2^ eq. 11]. It should be noted that to remove 
the impact of any shadowing processes, which are not accounted for in the K-/i fading model, 
the data sets were normalized to their respective local means prior to parameter estimation. To 
determine the window size for extraction of the local mean signal, the raw data was visually 
inspected and overlaid with the local mean signal for differing window sizes. For the D2D 
channel data, a smoothing window of 500 samples was considered. As we can quite clearly see, 
from Figs. [6ta) and (b), the envelope PDF of the K-/i fading model provides an excellent fit to 
the D2D data. To allow the reader to reproduce these plots, parameter estimates for all three 
measurement scenarios are given in Table |II| 

Using the parameter estimates obtained from the field trials. Figs. [6tc) depicts the estimated 
probability of SPSC versus 7 m for selected values of 'Je for the measured D2D channel. In this 
instance, it can be seen that the estimates for k of the main and eavesdropper’s channels are 
comparable and also greater than 0, suggesting that a dominant component existed for both. We 
also observe that the parameter estimates for /im and yUg are both quite close to 1, suggesting that 
a single multipath cluster contributes to the signal received by both Bob (node B) and Eve (node 
E) and thus this fading scenario is quite close to the Rician case. Eor this fading environment, 
we note that Pq increases as the average SNR of the main channel increases. Eurthermore, for 
a fixed 7 m, the channel becomes more susceptible to eavesdropping as the average SNR of the 
eavesdropper’s channel increases. 

B. Body Area Network Scenario 

The second set of measurements considered on-body communications channels operating at 
5.8 GHz as found in body area networks. The experiment was performed in the same seminar 
room discussed above which was unoccupied except for the test subject on which the on-body 
nodes were placed. Eor this scenario, to maximize coupling across the body surface, the antennas 
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Fig. 5. Seminar room environment showing the position of nodes A, B and E for (a) D2D scenario (b) on-body scenario. 

were mounted normal to the torso of an adult male of height 1.83 m and mass 74 kg. Node 
A was positioned at the front eentral chest region at a height of 1.42 m while nodes B and E 
were placed on the rear of the test subject at the central waist region at a height of 1.15 m and 
the right back-pocket at a height of 0.92 m, respectively. The measurements considered the case 
when the hypothetical BAN user walked along a straight line within the large room, covering a 
total distance of 9 m as shown in Fig. Ob). For the BAN scenario, a total of 19260 samples of 
the received signal power were obtained and used for parameter estimation. 

Figs.Od) and (e) show the empirical PDF of the signal envelope for Bob and Eve, again com¬ 
pared to the K-jj, PDF given in lf26l eq. 11]. Identical to the analysis of the D2D measurements, 
the optimum window size was determined from the raw channel data. In this case a smoothing 
window of 100 samples was used. Again, the K-jj, PDF was found to provide an excellent fit 
to the empirical data for both Bob and Eve. Interestingly for the BAN configuration considered 
here, the estimated n, parameter of the eavesdropper’s channel was greater than that of the main 
channel, while for the estimated /i parameters, the converse situation was true (Table HH). Fig. 
[^f) shows the probability of SPSC versus 7 m with selected values of for the measured BAN 
channel. It can be quite clearly seen that for fixed the probability of SPSC improves as the 
legitimate channel experiences a higher average SNR. 



































TABLE II 

Parameter Estimates for the k — Fading Model Obtained erom the Field Measurements 
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Eading Channel 

Km 


tm 

Ke 

Ae 

te 

D2D 

1.07 

0.91 

1.22 

1.11 

0.92 

1.19 

On-Body 

2.92 

0.75 

1.17 

3.60 

0.67 

1.17 

V2V 

5.02 

0.70 

1.04 

7.17 

0.60 

1.03 


C. Vehicle-to-Vehicle Scenario 

The third set of measurements eonsidered vehiele-to-vehicle communieation channels operat¬ 
ing at 5.8 GHz. The experiments were conducted in a business district environment in the Titanic 
Quarter of Belfast, United Kingdom. As shown in Fig. |71 the area consisted of a straight road 
with a number of office buildings nearby. For this particular scenario, nodes A, B and E were 
placed on the center of the dash boards of three different vehicles; namely, a Vauxhall (Opel 
in Europe) Zafira SRi, a Vauxhall Astra SRi and a Hyundai Getz. The initial positions of the 
vehicles are shown in Eig. |71 The measurements began when the vehicles that contained nodes 
A and B started approaching one another at a speed of 30mph. During these measurements the 
vehicle containing node E remained parked (with the driver still seated inside) on the side of 
the road as indicated in Eig. |71 It should be noted that all of the channel measurements made 
in this scenario were performed during off-peak traffic hours and were subject to perturbations 
caused by the driver, movement of the nearby pedestrians and other vehicular traffic. A total of 
56579 samples of the received signal power were obtained and used for parameter estimation. 

Figs.Og) and (h) show the empirical PDE of the signal envelope for Bob and Eve compared to 
the n-jj, PDE. Similar to the analysis of the D2D and BAN measurements, the optimum window 
size was determined from the raw channel data. The local mean for the V2V measurements was 
calculated over 200 samples. Erom these figures we can see that the PDE of the k-jj, fading 
model provides a very good approximation to the V2V data. Erom Table HH it can be seen that 
the estimates of n for the main and the eavesdropper’s channels are greater than 0 whilst the 
estimated /i parameters are less than 1, suggesting that a dominant component exists and that 
these channels suffer less from multipath. We also observe that the estimated k parameter of the 
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BAN: Rx-Node B 



0 0.5 1 1.5 2 2.5 3 

Received signal relative to the local mean 



V2V: Rx-Node B 



Received signal relative to local mean 



Fig. 6 . Empirical envelope PDF of node B and E compared to the K-fj, PDF given in I26I eq. 11] and the probability of SPSC 
versus 7 m with selected values of 7 ® for D2D, BAN and V2V channel measurements, respectively. 



Fig. 7. Satellite view of measurement environment showing the position of nodes A, B and E for V2V scenario. The vehicle 
with node E remained parked on the side of the road and was oriented such that it faced directly towards node A. Arrows 
indicate the front of the car. 
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eavesdropper’s ehannel was greater than that of the legitimate ehannel, the estimated yU parameter 
of the main ehannel was only marginally greater than that of the eavesdropper’s ehannel. Fig. 
mi) depiets the probability of SPSC versus 7 m with seleeted values of for the measured V2V 
ehannel. As observed for the previous two seenarios, the probability of SPSC is improved when 

7m > lE- 

VII. Conclusion 

Novel analytieal and elosed-form expressions for the probability of SPSC and SOP^ of the 
reeently proposed fading model have been presented. Speeifieally, the analytieal expres¬ 
sions have been derived for positive, real, i.n.i.d. K-fi ehannel variables in the presenee of an 
eavesdropper. We have also arrived at an exaet elosed form expression for the probability of 
SPSC for integer values of hm and he- Based on these results we have provided a useful 
insight into the behavior of SPSC and SOP^ as a funetion of parameters {km, ^m, 7m} 
and {ke, /xe, 7e}- The analytieal and elosed form expressions have been validated through 
reduetion to known speeial oases. As the n-fj fading model is a very general statistioal model 
that inoludes many well-known distributions, the new equations derived in this paper will find 
use in eharaoterizing the seereey performanee of several different fading ehannels (see Table 
IJ). Moreover, the results presented here will also find immediate applieation in the analysis 
of outage probability in wireless systems affeoted by eo-ehannel interferenoe and baokground 
noise, and the oaloulation of outage probability in interferenoe-limited seenarios. Finally, we 
have illustrated the utility of the new formulations by investigating the probability of SPSC 
based on real ehannel measurements oonduoted for a diverse range of wireless applioations suoh 
as eellular devioe-to-deviee, vehiele-to-vehiele and body eentrie fading ehannels. It is also worth 
highlighting that all of the expressions presented in this paper ean be easily evaluated using 
funetions available in mathematieal software paekages sueh as Mathematiea and Matlab. 
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Appendix A 
Proof of Equation (fT4l) 

From (fTSl) we have, 

Pout (Rs) = Pr [jM < (1 -I- '^e) — l] 

■e«S(l+7B)-l 


= / fiE (7e) 


/ f'lM (Tm) d'jM 


djE 


LXJ 

= I fiE ( 7 b) [PiM (1 + 'Je) - l)] d-fE 


Substituting dH) and dH) in (l26l) we obtain (fT^ . 


(26) 


Appendix B 

Proof of Proposition d] 

An analytieal expression for (fTVl) ean be derived by expressing the generalized Mareum Q- 
function and the modified Bessel funetion of the first kind aeeording to [1^ eq. 16] and 041 
as, follows. 


Qm P) 'y ^ 


(m + /, f) 


^ /!r(m + /)2'e^ 

v-\-2k 


tv(x) = y] 


f) 


fc =0 


klT {v + k + 1) 


(27) 


(28) 


where r(-) is the gamma funetion, r(-, •) is the ineomplete gamma funetion. By following the 
procedure proposed in O^ Lemma 1], the following upper bound is obtained for the truncation 
error of the Marcum Q- function representation in 

-I i-i (-i)'fefc() - k)i_^ [l - (-l)^e2“^ 


t ^ 


m-2 1-1 

1=1 k=0 


k\y/TT2 




+ Q {a+ b) + Q {b - a) - ^ ^ ^ 


2 0 


e 2 


(m + /) 2*e 2 
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By substituting (1271) and (l28l) in (fTTI) we obtain, 


SOP^ = 


Mg + l 

(3e ^ 

e"® + ^) 0 


CXD OO 


OO 


ME + l 

Pe ^ 

ctE 2 e"® 


X 


EE 

fc=0 i=0 


Hr()^r^«r(^„+T 2 l“" 1 7E"+‘-'e-'’'«r (f,„ + /, /3„e«»7E) <i7E 


where a = 1 /^m, b = 1/^e, Pm = (km + l)anM, Pe = + l)bfiE, c^m = ^M/iju and 

dE = ke/j^e- Notably, the integrals in (l29l) are identieal to |l3^ eq. 3.381,4] and |l32l eq. 6.455] 
given by. 


OO 


I X 
0 


V—1 


e 



and 


[Re jj > 0, Re V > 0] (30) 


OO 

/ {v, ax) dx 

0 


a"V (/i + v) 
fi{a + /3)^+" 


[Re {a + P) > 0, Re jj, > 0, Re {jj. + v) > 0] 


(31) 


We then substitute (l30l) and (l3TI) in (l29l) to obtain (IT^ . 


From (fT9l) we have. 


Appendix C 

Proof of Equation (I20l) 


Pq — Pr ( 7 m > 1e) 


= / fiE (7 b ) 


0 

OO 


/ ( 7 m ) d'jM 


SfE 


d'^E 


- / fiE ( 7 b ) [1 - -^ 7 m ( 7 b )] c ? 7 b - 
0 


(32) 


Substituting (l8|) and (Hj) in (l32l) we obtain (l20l) . 
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Appendix D 

Proof of Proposition [2] 

An analytical expression for the probability of SPSC is obtained by substituting (ITT!) and (1281) 
in (l20l) as follows, 

— l + 2fc /„ -.1 


<^£ + 1 oo oo 


— — k\T + k) l\T (/tm + 0 2'e"^ 

CXD 

X I (yUM + y PmIe) d'jE 

0 


(33) 


Notably, the above integral is identieal to 11321 eq. 6.455] given by. 


J{v,ax)dx = ^ ^^ 1^+v 2^1 + + 

0 ix{a + pf V a + p 


[Re {a + (3) > 0, Re p > 0, Re {p + v) > 0] (34) 


We then substitute (l34l) in (1331) to obtain (l2Tl) . 


Appendix E 

Proof of Proposition [3] 


From 1351 eq- 2.5], 


(A, 5; r) = f x^+^e-^—Jl^ (Ax) dx f y^+^e 


y^+B^ 


R (By) dy. (35) 


Using the definition of the Mareum (^-function given in (U) in (1351) . we obtain. 


Pf,^y{A,B]r) = l-A^fx'^~^^e ( = ^ R{Ax) Qy+i {B,rx) dx. 

0 


( 36 ) 
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Now letting x = and performing the neeessary transformation of variables we 

obtain, 


0 


-a //Te (1 + Ke) 


lE 


^1+1 


-Mg(l + W£;)7g 


7^2 6 T'-B 


X 4 Mi 


'2/iE (1 + i^e) 1e 


lE 


Qv+1 ( B, T ^ 


2/iE (1 + ke ) 7e 
1e 


diE. (37) 


Comparing (iJTl) with (l20l) and with the appropriate variable substitutions (see Proposition [3]), 
we obtain: 


Me+1 


/iE(l + Ke) 2 7 tMzl -^Ei^^EhE 

P^^^{A,B]r)=l --/7e 2 e 

2 7g 2 0 




= 1 - Pn. 


From [I35l eq. 3.16], we have. 


(38) 


P^,„ {A, B] r) = Po,o {A, B] r) + exp 


A^r + P"^r 


2.„-l 




m=—fi 


Br J 


R 


V + k 


X < ^ I I r“ 

k=i \ k + m 


v—k-\-lj^—v—k— 


j=i \ m 


(39) 


From 051 eq. 3.5], we have, 

Po,o {A, B]r) = Q 


B 


a/ 1 + r^’ vTTr^ 


(l + r^) ^ exp 


AV + P2 

’ 2(1 + r2) 


ABr 

1 _|_ j,2 


(40) 


Letting P = Pq^q (^4, P;r) and eombining (|38]) . (I39j) and (l40l) we obtain (l22l) . Note that ||35]| 


uses lower-ease symbols (a, b) and we use upper-ease symbols (A, B). This is beeause we define 


a = A- and h 


IM 


1 

lE ‘ 
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